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We present an analytical calculation of the spin- wave spectrum of the Jahn- Teller system LaTi03 . 
The calculation includes all superexchange couplings between nearest-neighbor Ti ions allowed by the 
space-group symmetries: The isotropic Heisenberg couplings and the antisymmetric (Dzyaloshinskii- 
Moriya) and symmetric anisotropics. The calculated spin-wave dispersion has four branches, two 
nearly degenerate branches with small zone-center gaps and two practically indistinguishable high- 
energy branches having large zone-center gaps. The two lower-energy modes are found to be in 
satisfying agreement with neutron-scattering experiments. In particular, the experimentally de- 
tected approximate isotropy in the Brillouin zone and the small zone-center gap are well reproduced 
by the calculations. The higher-energy branches have not been detected yet by neutron scattering 
but their zone-center gaps are in satisfying agreement with recent Raman data. 



PACS numbers: 71.10.- 



71.27.+a, 75.30.Ds, 75.50.Ee 



I. INTRODUCTION 



The orthorhombic perovskite LaTiOs has long been 
considered as a typical antiferromagnetic Mott insula- 
tor (Tn = 146 K). Albeit its rather small ordered mag- 
netic moment, 0.46— 0.57/is, 1,2 experimentally it seems 
not very different from a conventional Heisenberg anti- 
ferromagnetic insulator. Indeed, the spin-wave spectrum 
measured by Keimer et al. 1 is well described by a nearest- 
neighbor superexchange coupling having the value 15.5 
meV, accompanied by a weak ferromagnetic moment. 
The latter has been attributed to a small Dzyaloshinskii- 
Moriya interaction, of about 1.1 meV. The experiment 
reported in Ref. 2 has found that the antiferromagnetic 
order of LaTiC>3 has a G-type structure along the crys- 
tallographic a direction, while the ferromagnetic moment 
is along the c direction. 

Because of the unusually small ordered moment, it has 
been proposed 3 that perhaps the cubic Kugel-Khomskii 
Hamiltonian 4 could be taken as a starting point for a 
successful interpretation of LaTiC>3. However, this cubic 
model has some very unusual symmetries which inhibit 
the appearance of long-range magnetic order at non-zero 
temperatures. 5,6 At strictly cubic symmetry, the fivefold 
degenerate d-levels on the Ti ions are split by the crys- 
tal field of the oxygen octahedra into the lower threefold 
degenerate t 2g levels, (occupied in Ti by a single elec- 
tron) and the higher twofold degenerate e g levels. In real 
materials, those degeneracies are frequently lifted by the 
Jahn- Teller distortion. 

Figure 1 portrays the crystal structure of LaTiC>3 (the 
enumeration we use for the Ti sites is marked in the fig- 
ure). The unit cell contains four Ti ions, and the crystal 
has the symmetry of the space group Pbnm. The Jahn- 
Teller effect in LaTiC>3 is caused by the twisting of the 
Ti-0 bonds with respect to each other, i.e., by differ- 



ences between the 0-0 bond lengths which amounts to 
a deviation of certain O-Ti-0 bond angles away from 
90°. The distortion leads to a crystal field that splits 
the levels, 2 yielding a crystal-field gap of about 0.24 eV 
between the orbitally non-degenerate ground state and 
the first excited level, a value which has been confirmed 
by a study of photo-electron spectroscopy. 7 A compari- 
son of the optical conductivity and of Raman data shows 
that the lowest orbital excitation is centered at about 
0.25 eV. 8 This value is in excellent agreement with the 
estimate of the crystal-field splitting according to Ref. 
2. Furthermore, the non-degenerate ground-state orbital 
due to the crystal-field calculations given in Ref. 2 is 
consistent with the orbital order found in NMR measure- 
ments of the Ti-3d quadrupole moment. 9 The presence of 
orbital order at low temperatures has also been inferred 
from measurements of the dielectric properties and the 
dynamical conductivity. 10 

An explanation of the magnetism of LaTiOs, which is 
based on the crystal-field calculation given in Ref. 2, is 
presented in Ref. 11. The calculation included spin-orbit 
interaction on the Ti ions as well, and found accordingly 
that the superexchange coupling between neighboring Ti 
ions consists of the isotropic Heisenberg exchange, and 
the antisymmetric (Dzyaloshinskii-Moriya) and symmet- 
ric anisotropies, which appear as a result of the spin-orbit 
interaction. These anisotropies conspire together with 
the isotropic coupling to determine the magnetic order 
at low temperatures, shown in Fig. 2. By minimizing 
the magnetic energy of the classical ground state it was 
found 11 that the magnetic order of LaTiOs is primarily 
that of a G-type antiferromagnet, with the ordered mo- 
ment along the crystallographic a axis, accompanied by 
a weak ferromagnetic moment along the c axis, in good 
agreement with experiment. In addition, it was found 
that there is a small A-type moment of the spin compo- 
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nents along the b axis, which (although not yet detected 
in experiment) is allowed by the symmetry of the system. 




FIG. 1. The crystallographic structure of LaTi03. The 
ten Ti ions, which constitute the twelve inequivalent near- 
est-neighbor Ti— Ti bonds are enumerated. For simplicity, 
only the oxygen octahedra around the four crystallographi- 
cally inequivalent Ti sites are shown. La ions from two layers 
are depicted as small spheres. We use orthorhombic coordi- 
nates, in which the x, y, z axes are oriented along the crystal- 
lographic a, b, c directions. 

In this paper we calculate the spin-wave spectrum of 
LaTi03, which evolves from the magnetic ground state 
found in Ref. 11. Since the magnetic unit cell con- 
tains four sublattice magnetizations, the spin wave dis- 
persion consists of four branches. In the zero spin-orbit 
coupling limit, these four branches collapse into two 
branches, an acoustic mode and an optical one, which 
both are two-fold degenerate. Accordingly, we term 
the two branches which evolve from the (zero spin-orbit 
coupling) acoustic waves as 'acoustic modes', and those 
which evolve from the optical ones as 'optical modes'. At 
the Brillouin zone center, the energies of the two acous- 
tic branches do not vanish but have gaps, of magnitudes 
2.7 meV and 3.0 me V . These values are quite close to 
the zone-center gap of about 3.3 meV deduced from neu- 
tron scattering. 1 Furthermore, these two modes are ap- 
proximately isotropic in the Brillouin zone, again in good 
agreement with the neutron scattering experiment. 1 We 
find that the two optical modes are quasi-degenerate, 
having a zone-center gap of about 43.3 meV. These modes 
have not been detected yet by neutron scattering but are 



in good agreement with Raman data 8 where at low tem- 
peratures an excitation peak is seen, which is centered at 
about 40meV and which disappears at Tn- 

Our calculation employs linear spin-wave theory, 
which expresses the deviations of the spins from their 
ground state configuration in terms of Holstein-Primakoff 
bosons. Wc therefore begin our analysis by outlining in 
Sec. II the determination of that ground state configu- 
ration. We then continue to derive in Sec. Ill the spin- 
wave Hamiltonian, and to obtain the spin-wave disper- 
sion. Section IV contains a numerical study of the dis- 
persion curves, together with a detailed comparison with 
experiment. The summary of our results is presented in 
Sec. V. 




FIG. 2. The magnetic order of the Ti ions in the classical 
ground state of the effective spin Hamiltonian of the lattice. 
The ions are enumerated according to the sublattice to which 
they belong. 



II. THE MAGNETIC GROUND STATE 

The analysis of the magnetic structure of LaTiC>3, car- 
ried out in Ref. 11, involves several steps. First, a micro- 
scopic Hamiltonian containing the relevant interactions 
on the Ti ions and between nearest-neighbor Ti ions is 
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derived. Treating this Hamiltonian in perturbation the- 
ory, one then derives the superexchange interactions be- 
tween nearest-neighbor pairs of spins of the electrons in 
the ground-state orbitals. This effective spin Hamilto- 
nian is summed over the entire Ti-lattice, to obtain the 
magnetic Hamiltonian. Finally, one minimizes the result- 
ing magnetic energy of the system, to obtain the classical 
magnetic ground state. In this section we briefly review 
these steps. 

The derivation of the microscopic Hamiltonian starts 
from a point-charge summation of the static crystal field 
for the Ti ions, employing a full Madelung sum over the 
crystal. This determines the eigenenergies and the eigen- 
states of the static crystal field acting on each Ti ion, i.e., 
the crystal-field d states. The effective hopping between 
the d orbitals of nearest-neighbor ions via the interven- 
ing oxygens is then written in terms of a Slater-Koster 
parametrization of the Ti-0 hopping. The other inter- 
actions included in the microscopic Hamiltonian are the 
on-site Coulomb interaction and the on-site spin-orbit 
coupling on the Ti ions. In this way, the microscopic 
Hamiltonian pertaining to a pair of nearest-neighbor Ti 
ions (denoted m and n) takes the form 

Hmn = H mn + \mn- (1) 

Here 

TlO _ TTCf _|_ TTC /n\ 

n van n mn ' mw \ ) 

where H cf is the static crystal-field Hamiltonian, and H c 
describes the intra-ionic Coulomb correlations of a dou- 
bly occupied <i-shell. Because of the rather low symmetry 
of the system, in treating the Ti 2+ ions which appear as 
intermediate states of the exchange processes it is neces- 
sary to take into account the full on-site Coulomb inter- 
action matrix. 11 The other part of the Hamiltonian (1) 
is 

V mn = H Ifm + H mn , (3) 

in which H tun is the kinetic energy, described in terms 
of the effective hopping matrix, and H so is the spin-orbit 
interaction. This part is treated in perturbation theory, 
in order to obtain from the Hamiltonian (1) an effective 
spin Hamiltonian, pertaining to the spins of the two Ti 
ions, which acts within the Hilbert space of the fourfold 
degenerate ground state of the unperturbed Hamiltonian 
H°. 

The detailed perturbation theory presented in Ref. 11, 
is carried out to second order in H tun and up to second 
order in the spin-orbit coupling (scaled by the coupling 
strength A). This procedure yields a rich superexchange 
coupling between the spins of the non-degenerate crystal- 
field ground states of the Ti 3+ ions. For a pair of two 
nearest-neighbor Ti ions, the effective single-bond spin 
Hamiltonian is found to be 

h = 7 S -S 4-T) -fS x S + S • 4 s S 

lv mn ° ran m a ~ van y^m nj ' m 'ran u nJ 

(4) 



where J mn is the isotropic Heisenberg coupling (second- 
order in the tunnelling amplitudes, and independent of 
A), T> mn is the Moriya vector (second-order in the tun- 
nelling amplitudes, and first order in A), and A s mn is the 
symmetric anisotropy tensor (second-order in the tun- 
nelling amplitudes and in A) . As can be seen from Fig. 1 , 
there are 12 inequivalent nearest-neighbor Ti-Ti bonds 
in the unit cell of LaTiC>3. By the symmetry operations 
of the space-group Pbnm the magnetic couplings of all 
8 intra-plane bonds can be expressed in terms of those 
pertaining to the (12)-bond, and all 4 inter-plane ones in 
terms of those of the bond (13). 11 We list the numeri- 
cal values of the couplings in Table I. 12 (The Tables are 
given on page 12.) 

The magnetic Hamiltonian is found from the single- 
bond spin Hamiltonian (4), by summing over the en- 
tire Ti lattice. To this end, one decomposes the lattice 
into four sublattices, corresponding to the four inequiv- 
alent Ti sites of the unit cell (see Fig. 1). Although all 
four sublattice magnetizations are of equal magnitudes, 
their directions are all different. Denoting the sublattice 
magnetization per site by M i7 the macroscopic magnetic 
Hamiltonian is found to be 11 

//>.: = feM r M J +D|].(M l xM J )+M 1 .r tJ .M j ], 

ij 

(5) 

where ij runs over the sublattice pairs (12), (13), (24), 
and (34) of Fig. 1. Here, are the macroscopic isotropic 
couplings, Dfj are the Dzyaloshinskii vectors (to leading 
order in the spin-orbit coupling), which are the macro- 
scopic antisymmetric anisotropics, and are the macro- 
scopic symmetric anisotropy tensors (of second order in 
the spin-orbit coupling). The relations between those 
macroscopic couplings and the microscopic single-bond 
couplings are listed in Table II, and the inter-relations 
between the macroscopic magnetic couplings of different 
bonds, which are dictated by the space group symme- 
tries, are found in Table III. 

The minimization of the magnetic Hamiltonian (5) 
yields the magnetic structure shown in Fig. 2. Table IV 
lists the details of this structure, in terms of the canting 
angles ip and i9 according to Ref. 11. This structure is 
going to be the basis for the spin-wave expansion carried 
out in the next section. 



III. THE SPIN- WAVE HAMILTONIAN 

The deviations of the spins away from their directions 
in the classical ground state may be described in terms 
of Holstein-Primakoff boson operators. In our case, the 
system consists of four sublattices, which implies the in- 
troduction of four different bosonic fields, and, in turn, 
four branches in the spin- wave dispersion. 
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The first step in the standard calculation of spin-wave 
dispersions is the rotation of the local coordinates at each 
sublattice, i, such that the new z axis will point in the di- 
rection of the corresponding sublattice ground-state mag- 
netization, Mj. This rotation still leaves the freedom to 
choose the new local x and y axes, i.e., to rotate the new 
coordinate system around its z axis. Denoting the new 
local coordinate system by x' i7 y\ and z\ (i = 1,2,3,4), 
we find that the convenient choice for our purposes (ex- 
plained in Appendix A) is 



Ec' . a> .a' 
^m ^mn ^ni 

(mn) 



(9) 



where the primes denote the rotated quantities, 

S' = TT -S A' —TT -A -TT* (10) 



Zi = 



M ' 
M = I M,-| 



M, x 



Vi = 



mi 



%i ~ Vi x Zi, 



m. 



= V /(M?)2 + (M?)2. 



(6) 



Consequently, there is a local rotation matrix t/j, per- 
taining to each of the four sublattices, 

' m t -M?M*/mi -M?M?/m,i~ 
U t = Mf/m t -MV/rrn , (7) 
Mf Mf M? 

which rotates the orthorhombic into the local coordinate 
system. We now apply this local rotation to the spin 
Hamiltonian (4), re- writing it for convenience in short- 
hand notation 

Sm * A n 



We next introduce the Holstein-Primakoff boson 
fields 13 for each of the four sublattices. Since we consider 
only the Ti ions, it is convenient to use a coordinate sys- 
tem in which the Ti ions occupy the sites of a simple cubic 
lattice, of unit lattice constant (this picture is the ap- 
propriate one for comparing with the experimental spin- 
wave data, 1 as discussed in the next section). It is also 
convenient to use a coordinate system in which nearest- 
neighbor Ti ions are located along the axes (namely, to 
rotate the orthorhombic coordinates by —45° around the 
z axis, see Fig. 2). Denoting the boson fields of sublat- 
tice 1, 2, 3, and 4 by a R , 6 R , c R , and <i R , respectively, 
where R is the radius vector to Ti No. 1 in Fig. 1, the 
spin-wave Hamiltonian, in the harmonic approximation 
takes the form 



^mn ' S n > 



(8) 



v mn — / 
(mn) (mn) 

where A mn is the 3x3 superexchange matrix, compris- 
ing all three types of magnetic couplings. In the rotated 
coordinate system the spin Hamiltonian takes the form 



hsw = h\ 2 + ft-U + h\ 3 + h s 2i . 



(11) 



The first term here refers to the lower plane of the unit 
cell, the second to the upper plane, and the last two terms 
refer to the inter-plane couplings. Explicitly, 



h sl 



h sl 



13 



'24 



]T'{ [C 12 (l) + C 16 (l)](a R a R + 6 R 6 R ) 
R +a R [C 12 (2)(6 

R + fo R-2nJ + Cl6(2)(&R_ nx+nH + &R_„ x _ n J] 

+ a R [C 12 (3)(& R + 4_ 2n J +C 16 (3){bl l _ nx+ay +&k-n.-n,)]} 
]T'{ [C 34 (l)+C 38 (l)]( CR c R + 4d R ) 

R + 4[C34(2)(d R + d R _ 2n J +C 38 (2)(d R _ ax+ny +d R _ ax _ ny )] 

+ 4[cs4(3)(4 + 4_ 2n j + a 58 (3)(4-n, + n„ + 4-n.-=,)]} + h - c - 

{ c i3(1)(4«r + c r c r) 
R + C 13 (2)4(c R + c R _ 2n J + C 13 (3)a R (c R + 4_ 2n J } + h. c, 

c 24 (i)(46 R + 4d R ) 

R + 0-24(2)4^ + d R _ 2n J+ C M (3)4(4 + 4_ 2n J}+ h.c. 



(12) 



Here n Q is a unit vector along the a-direction, where a — x, y, z, and the Ti ion marked by 1 in Fig. 1 is at the origin. 
The summations then extend only over the Ti ions No. 1 in each unit cell (this is indicated by the prime on the 
summation symbols). The coupling coefficients C mn {€) in Eq. (12) are given by combinations of the superexchange 
matrix elements (A' mn ) al3 , 

Cmn(l) — — ftAnn) ' 

Cmn(2) = \ \{^ mn ) xx + {A' mn ) yv + i((A' mn ) yx — {A' mn ) Xy ^ , 

C„ m (3) = \ [{A' mn r - (A' mn ) yv + *((A' mn ) yx + (A' m n) xv )} ■ (13) 
In writing down Eq. (12), we have omitted constant terms. The transformation to the Holstein-Primakoff operators 
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yields also terms which are linear in the boson fields; these vanish upon summing over all bonds (see Appendix A) 
due to the proper choice of the local coordinate system. 

Our magnetic unit cell is spanned by the vectors (1,1,0), (1,-1,0), and (0,0,2), and the corresponding magnetic 
Brillouin zone (MBZ) is defined by 



\q x + q y \ < 7T, \q z \ < -. 



By introducing the Fourier transforms of the operators, 



n X] 

qGMBZ 



J R 



J_ "V P *q(R-+n*)!,t 
n 2-^i q' 



qGMBZ 



, J_ jt _ / J. V p »q(R+n*+n*)jt 

R — \/ n Z_^i e L q' U R — y N Z_^i e u q' 

qeMBZ qeMBZ 



where N is the total number of magnetic unit cells, the spin-wave Hamiltonian (11) becomes 

hsw = ^2 /isw(q), 



(14) 



(15) 



(16) 



where 

/isw(q) = Ci (a q a q + 6 q 6 q + c q c q + <i q d q ) 

C| (cos q x + cos g y ) (a q 6 q + c q d q ^ + h. c. + C^ cos g z (a q c q + 6 q d q ) + h. c. 



+ 
+ 



(Cj cos g x + C|* cos 9y ) (a q 6^ q + c q dL q 



h. c. 



Cf cos g z a q cL q + Cg 1 * cos ^6 q rf f _ q + h. cj . (17) 



The coefficients appearing in this equation are linear combinations of the previous coefficients C mn (i) (see Appendix 
A), 

C 1 = 2C 13 (l)+4C 12 (l) = Cr, Ci = 2C 13 (2) = Ct , C\ = 2C 12 (2), C3 1 = 2C 13 (3), cj = 2C 12 (3). (18) 

These are related to the original spin-coupling coefficients of Eq. (4), but arc not reproduced here explicitly, since 
their expressions are very long. 

The spin-wave dispersion pertaining to the Hamiltonian (16) is calculated in Appendix B, leading to the result 

n?(q) - (C\ - C£ cos q z f - \Cf\ 2 cos 2 q z + \c\ | 2 (cos fe + cos^) 2 

- |C| cosq x + Cj* cos^| 2 - (cos,q x + cos q y ) W (cos q z ) , 
ft 2 (q) = ft 2 (q + Q), with Q = (0,0, tt), 



where 



2 (q)=0 2 (q+Q'), with Q' = (tt, tt, 0), 

ft 2 (q) = f7 2 (q + Q"), with Q" = Q + Q' = (^,tt,^), 



iy 2 (cos^) =4[(d - C^cosg,) 2 - l^pcos 2 ^] [|C|| 2 - ( C3+ 2 C3 * )' 



(19) 



+ 



(C^C, 1 + CtCf) cosq z + (d - cos^)^ 1 + Cf ) 



(20) 



Each of the branches has tetragonal symmetry, i.e., 
tti(q x ,q y ,q z ) = ^i{q y ,q x ,qz) = V t (-q x , q y , q z ) = 
tti{qx,-q y , qz) = ^(q x ,q y , -q z ). 

Equations (19) contain our final result for the spin- 
wave spectrum of LaTiC>3. Evidently, the details of the 



spectrum can be obtained only numerically: One has to 
write the spin-wave coefficients, Eqs. (18), in terms of 
those appearing in Eqs. (13), and express the latter via 
Eqs. (8) and (10) in terms of the original coefficients of 
the spin Hamiltonian (4) using the values listed in Table 
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I. These results are then used in constructing the dis- 
persion. We carry out this procedure in the next section, 
confining ourselves to the wave vectors explored in the 
neutron scattering and Raman experiments, respectively. 

When the spin-orbit coupling A is set to zero, the co- 
efficients appearing in Eqs. (19) simplify to 



C\ — 2Ji2 - 

-■^13, 



13, 

cl 



c\ =0, 



'12, 



where Ji 2 is the isotropic in-plane Heisenberg coupling, 
and J 13 is the Heisenberg coupling between planes. In 
that case [note that cos q z = | cos q z | in the Brillouin zone, 
see Eq. (14)] 



n?(q) = nl(q) - (2J 12 + J 13 ) 2 

- (j 12 (cosq x + cosq y ) + J 13 \cosq z \Y 



(22) 



while the expression for J7| (q) = f2| (q) is obtained upon 
changing cos q x + cos q y to — cos q x — cos q y . At the zone 
center f2i and f2 2 vanish, while O3 and O4 have a gap 
equal to V8J12J13. Obviously, in the absence of spin- 
orbit coupling the magnetic unit cell includes only two 
sublattices (in that case, sublattice 1 and sublattice 4 
can be combined into one sublattice, and so can sublat- 
tice 2 and sublattice 3 in Fig. 2). The Brillouin zone 
corresponding to this smaller magnetic cell is twice as 
large as the one of Eq. (14). By "folding out" the optical 
mode into this larger Brillouin zone, one reproduces the 
usual gap-less dispersion of the pure Heisenberg model. 
At finite values of the spin-orbit coupling all modes have 
gaps at the zone center, but those of Oi and f2 2 are 
much smaller than the ones of the other two modes. For 
this reason, we term the f2i(q) and the f2 2 (q) branches 
'acoustic modes' and r^(q) and f2/i(q) are referred to as 
optical modes. Optical spin-wave modes have been de- 
tected, for instance, in bilayer cuprates. 14-16 



IV. NUMERICAL RESULTS FOR THE 
SPIN- WAVE DISPERSION 



For the model parameters we use, it turns out that the 
two acoustic branches as well as the two optical branches 
are nearly degenerate. The reason is the smallness of the 
angle ip, which leads to an additional translational sym- 
metry which is nearly fulfilled by the classical ground 
state. This "quasi" -symmetry corresponds to the trans- 
lation by the vector R44 which connects the Ti ions No. 1 
and 4 (see Figs. 1 and 2). For p=0 this symmetry is ex- 
act, and the magnetic unit cell contains only two ions. 
In that case the spin-wave dispersion consists of two 
branches. As we have a small deviation from this ideal 
case, we obtain two pairs of quasi-degenerate branches. 



A. Comparison of the acoustic branches with 
neutron scattering data 

We begin our discussion here by recalling the experi- 
mental results of Rcf . 1 . The authors of Rcf . 1 have fitted 
their neutron scattering data with an isotropic single- 
branch spectrum parametrized as 



(21) fi(q) 




-(cosq x + cosq y + cosq z ) 2 . (23) 



This assumes an isotropic Heisenberg coupling, J, for the 
entire Ti lattice, namely, the same coupling for the bond 
(12) and the bond (13) of Fig. 1, and introduces a zone- 
center spin-wave gap, A. The experimentally determined 
values of these parameters are 



J= 15.5 ±1.0 meV, A = 3.3 ± 0.3 mcV. 



(24) 



In the following, we compare the fitted function, Eq. 
(23), with the acoustic branches fii(q) and fi 2 (q). 

Although the symmetry of our spin- wave Hamiltonian 
allows for two acoustic modes, the resolution of the dis- 
persion measurements, which amounts to about 10 % at 
any given point q in the Brillouin zone, 17 is insufficient 
to resolve the two branches. To demonstrate this point, 
and to compare in detail the experimental findings with 
our expressions, we proceed as follows. Firstly, we aver- 
age the Heisenberg couplings pertaining to the different 
bonds (calculated in Ref. 11) over the six Ti-Ti bonds in 
which each Ti ion is participating, 



4 J 12 + 2 J; 



13 



G 



= 15.89 mcV. 



(25) 



Clearly this value agrees with the experimental one given 
in Eq. (24), within the accuracy of the experiment. Sec- 
ondly, we calculate the zone-center gaps as found from 
our calculation. Following the numerical procedure out- 
lined at the end of the previous section, we find 



Ai = fii(O) = 2.71 mcV, 
A 2 = O 2 (0) = 2.98 meV. 



(26) 



We have found that the splitting between the two calcu- 
lated acoustic branches reaches its maximum at the zone 
center, where 



=91.14%. 

A 2 



(27) 



This discrepancy is within the uncertainty of about 10 % 
of the measured spin- wave energies of Rcf. 1. 

Away from the zone center the two acoustic branches 
are quasi-degenerate. We estimate the tetragonal 
anisotropy of the dispersion by comparing the dispersions 
at wave vectors q = (n/2, 0, 0) and q = (0, 0, 7r/2), 



fii(f,0,0) 



91.34%, 



n 2 (o,o,j) 
n 2 (f,o,o) 



91.29%. (28) 



6 



This implies that the tetragonal anisotropy is also less 
than the uncertainty of the measured spin-wave energies. 
The calculated dispersions along selected directions in 
the Brillouin zone are depicted in Figs. 3, together with 
the optical branches which we will discuss in Sec. IV B 
and the experimental dispersion computed from Eq. (23) . 
The agreement between the acoustic branches and the 
experimental dispersion is satsifying. 

It is harder to infer the experimentally quoted value 1 of 
the Dzyaloshinskii-Moriya interaction, 1.1 meV, (which 
does not agree well with our values for the Dzyaloshin- 
skii vectors, see Table I), from the calculated disper- 
sion. We therefore attempt to estimate the effects of the 
two types of anisotropics, antisymmetric and symmetric, 
on the spin- wave dispersion by analyzing two cases: (i) 
Switching off all antisymmetric anisotropics, D mn = 0, 
(all other terms are accounted for according to their cal- 
culated values, see Table I) and (ii) switching off all sym- 
metric anisotropies, A s mn = 0, while keeping the con- 
tributions of the antisymmetric ones. In both cases we 
examine the spin canting, i.e., the ground-state config- 
uration of the magnetization, and the zone-center gap 
of the dispersion. (The dispersion away from the zone 
center is dominated by the Hciscnberg couplings.) 

(i) In the absence of the Dzyaloshinskii-Moriya inter- 
action, the canting practically disappears. We find that 
the canting angles almost vanish, 



-0.04°, = 0.00°, if D, 



0. 



(29) 



However, the zone-center gap is enhanced compared to 
its actual values, Eq. (26), 



Aj = A 2 =4.73meV, if D ? 



0. 



(30) 



(ii) In the absence of the symmetric anisotropies the 
spin canting is almost the same as given in Table IV, 



^ = 1.47°, tf = 0.80°, if A s mn = 0. 



(31) 



Switching off continuously the symmetric anisotropies, 
we find that the zone-center gap first closes and then even 
becomes imaginary as the symmetric anisotropies ap- 
proach zero. This unphysical result shows that one is not 
allowed to consider only the antisymmetric anisotropics 
resulting from the spin-orbit interaction, without includ- 
ing the symmetric ones as well. Indeed, as has been 
already pointed out in Refs. 18 and 19, a systematic 
treatment of the effect of the spin-orbit interaction on 
the spin couplings must include both anisotropies. They 
both contribute to the magnetic energy terms of the same 
order in the spin-orbit coupling parameter. 

Comparing these two fictitious cases, we conclude 
that the spin-canting is dominated by the antisymmetric 
anisotropies, while the zone-center gap of the dispersion 
is governed by the symmetric anisotropies. It is therefore 
a somewhat questionable procedure to deduce the anti- 
symmetric anisotropy of the spin coupling from the spin- 
wave dispersion, taking into account only the Moriya vec- 
tors, as has been done in Rcf. 1. This is again related to 



the fact that both the Dzyaloshinskii-Moriya interactions 
and the symmetric anisotropies induced by the spin-orbit 
coupling appear in the same order in the magnetic energy 
and in the spin- wave dispersion. 20 

The manner by which the various anisotropic spin cou- 
plings, in a low-symmetry system like LaTiOs, can be de- 
duced from an experimentally obtained spin-wave spec- 
trum therefore remains unsettled. In our case, the spin- 
wave Hamiltonian, Eq. (17), depends on 8 parameters 
[note that some of the coefficients, Eqs. (18), are com- 
plex]. Furthermore, even the knowledge of these 8 param- 
eters does not suffice in our case to trace backwards the 
parameters of the spin Hamiltonian, Eq. (4). The reason 
being that the coefficients involving the matrix elements 
(A' mn ) xz and (A' mn ) yz [see Eq. (13)] disappear altogether 
from the spin- wave Hamiltonian (see Appendix A). The 
conclusion is that it is possible to use certain numeri- 
cal values for the various types of spin couplings and to 
investigate their consistency with the experimentally de- 
tected spin- wave dispersion (as done above). However, an 
unequivocal deduction of spin-coupling parameters from 
spin-wave spectra is not possible due to the low symme- 
try of this system. 



B. The optical branches 

The two calculated optical branches, depicted in 
Figs. 3, are practically indistinguishable in the entire 
Brillouin zone. Their zone-center gaps are 



A 3 = fi 3 (0) = 43.32 meV, 
A 4 = Ot(0) = 43.34 meV. 



(32) 



So far, branches with such a large zone-center gap have 
not been detected by neutron scattering. 1 Possible rea- 
sons are: (i) The signal in the energy range of the op- 
tical branches has a rather low intensity (as compared 
to the lower energy regions); (ii) The spin- wave signal 
in this energy range is accompanied, and possibly is hid- 
den, by phonon excitations. 1 However, despite of these 
two problems, in principle it might be possible to de- 
tect the dispersion of the optical branches by neutron 
scattering. 17 Our prediction is that the dispersion of the 
optical modes will be qualitatively different from that of 
the acoustic ones. These modes will not have the ap- 
proximate isotropy of the acoustic modes, but will show 
a larger tetragonal anisotropy. We find 



n 3 (o,o,s) M Ot(o,o,f) M , 

2 - 70.47 %, - i l[ = 70.44 %. (33) 



o 3 (f,o,o) • ' n 4 (f,o,o) 

These relations can serve as a further check of our model. 
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(a) q= f(l,l,l)^ 



(b)q = 7r(l,0,0)£ 
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FIG. 3. The spin-wave dispersion along selected directions 
in the magnetic Brillouin zone. We use pseudo-cubic coordi- 
nates, in which the Ti ions No. 1 and 2 are located along the 
x axis. Panels (a)-(c) show (in black) the four branches fi,(q) 
of the calculated dispersion and (in grey) the single branch 
n(q) which has been fitted onto neutron scattering exper- 
iments, Eq. (23). The acoustic branches fli(q) and ^(q) 
are quasi-degenerate, such that away from the zone center no 
splitting between them can be seen. The optical branches 

(q) and fLt(q) are practically indistinguishable over the 
entire Brillouin zone, (a) The dispersion along (1,1,1). This 
direction is chosen because the experimental paper on the 
neutron scattering contains a plot along this direction where 
the measured points of the dispersion are shown. 1 Though 
the calculated acoustic branches of the dispersion give slightly 
lower energies at the zone-center than the fitted function and 
slightly higher energies at the zone edge, these deviations are 
within the uncertainty of the measurement and hence, the 
agreement of our calculated dispersion with the the measured 
points and with the fitted function is satisfying. The splitting 
of the calculated acoustic branches at the zone center is too 
small to be resolved in the experiment. From panels (b) and 
(c) one can see that the tetragonal anisotropy of the calculated 
acoustic branches is rather small. The agreement between the 
acoustic branches and the neutron scattering data is satisfy- 
ing also along the (1,0,0) [panel (b)] and (0,0,1) [panel(c)] 
directions. 
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In contrast to the absence of experimental evidence 
for the optical modes in the neutron scattering exper- 
iment, Raman spectroscopy 8 at low temperatures docs 
show a pronounced peak centered at about 40meV. This 
energy is consistent with our calculated optical branches 
3 (q) and (q). In Raman spectroscopy only the zero 
wave vector excitation of the optical branches can be ob- 
served. In principle, Raman spectroscopy is only sensi- 
tive to S z —0 excitations but this selection rule can be 
broken by the spin-orbit coupling. The Raman peak dis- 
appears at the Neel temperature, giving evidence for a 
magnetic origin. Studying spin-wave energies in Raman 
spectroscopy might be subject to similar difficulties as 
neutron scattering when it comes to the phonons' role. 
Since the pronounced peak at about 40 meV has a very 
large intensity, its explanation may well have to include 
the coupling to lattice modes, in addition to the optical 
spin-wave modes. 
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(c) q= f(0,0, IX 



[meV] 




and are approximately isotropic over the Brillouin zone. 
The dispersion and the gap of these two modes are shown 
to reproduce the experimental data of the neutron scat- 
tering experiment carried out on LaTiOs. 1 The quasi- 
degenerate modes belonging to the second pair have a 
large zone-center gap, about 43meV, and their disper- 
sion shows sizeable tetragonal anisotropy in the Brillouin 
zone. While not yet detected in neutron scattering ex- 
periments, (perhaps for technical reasons as indicated 
above), the zone-center gap of these modes is consistent 
with Raman data. 8 

Our spin-wave dispersion is calculated on the basis 
of the detailed low-temperature magnetic structure of 
LaTi03, which we have analyzed in a previous paper. 11 
There, we have used the experimentally verified orbital 
ordering in this system, to develop the superexchange in- 
teraction between nearest- neighbor Ti ions. As detailed 
in Ref. 11, and summarized in the Introduction section 
above, the complicated magnetic structure that we have 
obtained, which involves a predominant G-typc antifer- 
romagnetic order along the a axis and a canted ferro- 
magnetic one along the c axis, agrees beautifully with 
all available experimental findings. In view of the good 
agreement we have found in the present study with the 
neutron and Raman scattering data, it might be con- 
cluded that our analysis has yielded a detailed under- 
standing of the magnetism in LaTiO,3. In addition, we 
have indicated above a rather detailed prediction regard- 
ing the behavior of the higher-energy modes. We hope 
that these will be studied experimentally, and will be 
compared with our calculations. 



V. SUMMARY 
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APPENDIX A: THE CHOICE OF THE 
ROTATION MATRIX 



As is mentioned in the text, the local coordinate sys- 
tem given in Eq. (6), in which the local z axis points 
along the direction of the moment in the classical ground 
state is still ambiguous in that it can be rotated arbitrar- 
ily around its z-axis. Here we show that our choice, Eqs. 
(6) and (7), leads to a considerable simplification in the 
calculation of the spin-wave Hamiltonian. 

Let us suppose that the local coordinate system of Eq. 
(6) at each of the four lattices sites is further rotated 
around its local z-axis by an angle pi, (i = 1,2,3,4). 
The rotation matrices Ui of Eq. (7) are then replaced by 



Ui 



ipi 



cos pi 
sin pj 




-sin pj 

COS Pi 





(Al) 



and the corresponding superexchange matrices, Eq. (10), 
are transformed accordingly as 



Al 



mnp m p n 



u. 



mp r , 



■ A' -TT 



where 



A' = A' 



(A2) 



(A3) 



is the superexchange matrix of Eq. (10). 

The arbitrary rotations described above will modify 
the coefficients C mn (£), Eqs. (13), appearing in the spin- 
wave Hamiltonian. Denoting these modified coefficients 
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by C mnPmPrl (e), such that C mn0 o(£) = C mn {l), we find in Table III) 
the following inter-relations (using the symmetries listed 

Ci 2 p lP2 (l) = c 34p3 p 4 (i) = c 16piP2 (i) - C 38p3P4 (l) = C 12 (l) = C? 2 (l), 

c 13P1P3 a) - c 24P2P4 (i) = c 13 (i) = c 13 (i), 

C 12[16]P1P2 (2) = C 12[16] (2)e < ^-" 1 ), C 34[38]p3P4 (2) = C 34[38] (2)e^-^ , 
^i2p lP2 (2) = C 16piP2 (2) = C 34p3P4 (2) = C 38p3P4 (2) if p 4 - P3 = P2 - Pi, 
C 13piP3 m = ^4 P2P4 (2) = C 13 (2) = CT 3 (2) 



C 



12[16]pip 2 l 



(3) = C 12[16] (3)e-^+^), C, 



34[38]p 3 P4 ' 



(3) = C 34[38] (3) e -^+.4) 



Ci2p lP2 (3) = C' 1 * 6piP2 (3) = C 34p3P4 (3) = C 38p3P4 (3) if Pl + p 2 = p 3 + p 4 , 
C 13pi p 3 (3) = CnWe-^+ri, C 24p2P4 (3) = C 24 (3) e -^+^, C 13 (3) = C 2 * 4 (3). 



(A4) 



It is thus seen that with the choice employed in Eq. (6), namely pi = (i = 1,2,3,4), the following relations are 
obtained 

Ci 2 (l) = Ci 6 (l) = C 34 (l) = C 38 (l), C 13 (l) = C 24 (l), 2C 12 (1) + C 13 (l) = 
C 12 (2) = Ci 6 (2) = C 34 (2) = C 38 (2) = C 13 (2) = C 24 (2) = hc£ = hot 



C 12 (3) = C 34 (3) = Cr 6 (3) - C* 8 (3) = \c\, C 13 (3) - C 2 * 4 (3) = V. 



(A5) 



Here we have introduced the coefficients C\, c\ ,A ~ , and that are used in our spin- wave Hamiltonian, Eq. (17). 

As is mentioned in the text, the Holstein-Primakoff transformation gives rise to terms linear in the boson operators. 
The coefficients of these terms are 



C mn (4) = -((A' mn ) xz +i(A' mn r)- 



(A6) 



When summed over all single-bond contributions, these coefficients vanish. For example, the absolute value of the 
coefficient of the boson operator is 



2|C 12 (4) + C 16 (4) + C 13 (4)| = y{A[r + A[* z + A'^} 2 + W + *J" + KT? = 0- 



(A7) 



Employing Eqs. (10), we have written each of the terms appearing in the square root explicitly, and verified that they 
both vanish. A similar argument prevails for the other coefficients of the linear terms. 



APPENDIX B: THE SPIN- WAVE DISPERSION 



In order to obtain the spin- wave dispersion resulting from the Hamiltonian (17), it is convenient to first introduce 
a short-hand notation for this Hamiltonian. To this end we write 



^sw(q) = £(?^(q)^(q)C(q) + ^^(q)^( q )^t ( _ q) + lw£(q)^(q)£„(-q)), 



(Bl) 



where 



€(q) 



, € f (q)=[4> b l c l 4 



q J ' 



and the Hamiltonian matrices are conveniently written in the form 



n a ( q ) = 



Hi Ti.2 
Hi Hi 



, Hi 



1 C\(oasq x + cosq y ) 

c\ * (cos q x + cos q y ) Ci 



, H2 — C2 cos q z 



1 
1 



(B2) 



(B3) 
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and 



TI3 TI4 
H4 TL3 



H3 = (Cg cosq x + Cj* cosg y ) 



1 

1 



Hi = cosq z 



C3 
ct 



(B4) 



[Note that C\ and C 2 are real, see Appendix A.] 

Let us now denote the boson fields in which the Hamiltonian (17) is diagonalized by r^(q), i — 1,2,3,4. These 
fields are related to the original ones, &(q), by the general linear transformation 



r,(q) = J2 ^(q)^-(q) - E Q«(<0£j(-q)> 



(B5) 



with 



E ( P « (<0 P »j (<0 - 0« (<l)Onj (0)) = ^ E (<l)Qnj ("<0 + Q/j Wry ("O)) = 0, 



(B6) 



for the r fields to obey the boson commutation relations. In order that the r fields will represent normal modes, they 
have to satisfy 



T £ (q),/isw(q) =0 £ (q)^(q), 



(B7) 



where f^(q), £ = 1,2,3,4 are the eigenfrequencies of our spin-wave Hamiltonian. Inserting Eqs. (B5) into Eq. (B7), 
and equating the coefficients of £ and £' on both sides, we obtain 



W«(q) = E(^n(q)K,(q) + Q,„(q)<(q)), 

n 

-fi/(q)Q«(q) = E(^n(q)KKq) + ^JqX-(q))- 



(B8) 



Identifying . = v* as "vector number ^ whose entries are j" , and similarly for Qij = Uj we arrive at the equations 

n e v e = n a *v l + n b *u l , -n e u l = n a u l + wV, (B9) 

where we have dropped the explicit q dependence for brevity. From the first of Eqs. (B6), we have 

\v e \ 2 - \u e \ 2 = 1, (BIO) 

where u l and v l are 4-dimensional vectors. 



We split the 4-dimcnsional vectors u and v into two 2-dimensional vectors, u = (u\, u 2 



(v{,v 2 ), an d write 



explicitly Eqs. (B9), using the definitions (B3) and (B4). The resulting equations may be arranged in the form 

- u 2 ) = (H 1 - H 2 )(u 1 - u 2 ) + (H 3 - H A )(v 1 - v 2 ), 

- v 2 ) = (HI - n 2 )( Vl - v 2 ) + (H* 3 - H* 4 )( Ul - u 2 ), 
+ u 2 ) = {H x + H 2 )( Ul + u 2 ) + (H 3 + W 4 )(« 1 + v 2 ), 

+ « 2 ) = (WI + W 2 )(«i + « 2 ) + («5 + + "2), (BH) 

where we have also dropped the index I for brevity. It is thus seen that there are two types of solutions: Either 
u 1 = u 2 and v 1 — v 2 , in which case the first couple of equations is trivially satisfied, and it is needed to solve just the 



second pair of equations, or vice versa: u 1 = —u 2 and v x 



and then the first pair of equations has to be solved. 



However, the only difference between the first pair of equations and the second one are the signs appearing in front 
of H 2 and H4. Glancing at Eqs. (B3) and (B4) reveals that these signs are determined just by cosq z . Therefore, it 
suffices to solve one pair of equations, and the solution of the second is obtained by simply changing the sign of cosq z . 
Focusing on the first option, we find that two of the eigenfrequencies are determined by 



dct 



Ci + cos q z + n 

C 2 *(cosq x + cosq y ) 

Cg 1 * cos q z 
Cj* cosq x + cj cos<7j, 



C 2 (cosq x + cosq y ) 



C3 cos q z 



C 3 cos q x + Cg* cos q y 



C\ + C 2 cos q z + Q C 3 cos q x + C 3 cos q y 
C|* cos q x + Cj cos q y C\ + C 2 cos q z — f2 
Cg 1 cos q z C| (cos q x + cos q y ) 



cos q z 



C 2 *(cosq x + cosq y ) 
C\ + C 2 cos q z — il 



= 0. 



(B12) 
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The modes ^(q) and fLi(q) [see Eqs. (19)] are the 
positive roots of the fourth-order polynomial in f2 given 
by Eq. (B12). The other two eigenfrequencies are found 
by changing the sign of cos q z . 

TABLE I. The single-bond spin-exchange couplings 
(in meV). The symmetric anisotropics are given as 

Ail _ f A xx Ayy Azz \ j lod _ ( Ayz Axz Axy \ r 

the diagonal and off-diagonal entries, respectively. 

Heisenberg couplings 
J12 = 17.094, Ji 3 = 13.484 
iva vectors 



TABLE III. Symmetries of the magnetic Hamiltonian due 
to the space group. The relations for the anisotropic couplings 
are abbreviated as follows: (+,+, +)i2 = (— ,+, +)ie means 
that D? 2 = {-Df 6 ,Df 6 ,Df 6 ), etc. Due to the glide planes, 
the Dzyaloshinskii vectors of the planar bonds have vanish- 
ing x components, and the respective symmetric anisotropies 
have vanishing xz and xy entries. Because of the mirror 
planes, the Dzyaloshinskii vectors of the inter-planar bonds 
have vanishing z components and the respective symmetric 
anisotropies have vanishing yz and xz entries. 

Isotropic couplings 

Il2 = I'M, Il3 = 124 



D12 = (2.260, -0.884, -0.893), Di 3 = (-2.207, 0.377, 0) 



Symmetric anisotropies 



A? 2 = (0.131,0,0), A? 3 = (-0.027,0,0), 
A?a = (0, -0.077, -0.061), Ali = (0,0, -0.052) 



Dzyaloshinskii vectors 



(0, +, +)l2 = (0, -, +) 3 4, ( + , +, 0)l3 = 



,0)24 



Macroscopic symmetric anisotropies 



(+, 0, 0)i2 = (-, 0, 0) 34 , (0, 0, +)i 3 = (0, 0, -) 24 



TABLE II. The macroscopic couplings of the sublattice 
magnetizations in terms of the microscopic single-bond spin 
couplings. For instance, J12 = J12 but /13 = J13/2, because 
the coordination number of a Ti ion is 4 in the planes and 2 
between the planes. 

Isotropic couplings 
I12 = J12, /13 = 
Dzyaloshinskii vectors 

D? 2 = (0,D y 12 ,D z 12 ), Dg ; = |Di3 

Macroscopic symmetric anisotropies 

rf 2 = a? 2 , ril = (A^,o,o), r 13 = ia 13 



TABLE IV. The structure of the magnetic order, charac- 
terized by the sublattice magnetizations Mi in the classical 
ground state (normalized to M), in terms of the canting an- 
gles (fi and We use orthorhombic coordinates, in which 
the x,y,z axes are oriented along the crystallographic a,b,c 
directions. 

x components: G-type 
-Ml = M% = M§ = -M% = M cos (p cos 1? 

y components: A-type 
-Ml = -Ml = Ml = M\=M sin ip cos ■& 
z components: ferromagnetic 
Mf = Mi =Mi = Mi = Msintf 
Calculated values of the canting angles 
ip = 1.42°, ■& = 0.80° 



1 B. Keimer, D. Casa, A. Ivanov, J. W. Lynn, M. v. Zim- 
mermann, J. P. Hill, D. Gibbs, Y. Taguchi, and Y. Tokura, 
Phys. Rev. Lett. 85, 3946 (2000). 

2 M. Cwik, T. Lorenz, J. Baier, R. Miiller, G. Andre, F. 
Bouree, F. Lichtenberg, A. Freimuth, R. Schmitz, E. Miiller- 
Hartmann, and M. Braden, Phys. Rev. B 68, 060401(R) 
(2003). 

3 G. Khaliullin and S. Maekawa, Phys. Rev. Lett. 85, 3950 
(2000). 

4 K. I. Kugel and D. I. Khomskii, Fiz. Tverd. Tela (Leningrad) 
17, 454 (1975) [Sov. Phys. Solid State 17, 285 (1975)]; Usp. 
Fiz. Nauk 136, 621(1982) [Sov. Phys. Usp. 25, 231 (1982)]. 

5 A.B. Harris, T. Yildirim, A. Aharony, O. Entin-Wohlman, 
and I. Ya. Korenblit, Phys. Rev. Lett. 91, 087206 (2003). 



6 A.B. Harris, T. Yildirim, A. Aharony, O. Entin-Wohlman, 
and I. Ya. Korenblit, Phys. Rev. B 69, 035107 (2004). 

7 M. Haverkort, Z. Hu, A. Tanaka, G. Ghiringhelli, H. 
Roth, M. Cwik, T. Lorenz, C. Schuessler-Langeheine, S. V. 
Streltsov, A. S. Mylnikova, V. I. Anisimov, C. de Nadai, 
N. B. Brookes, H. H. Hsieh, H.-J. Lin, C. T. Chen, T. Mi- 
zokawa, Y. Taguchi, Y. Tokura, D. I. Khomskii, and L. H. 
Tjeng, cond-mat/0405516 (2004). 

8 A. Gofiling, M. Griininger, B. Keimer, and C. Ulrich, private 
communication. 

9 T. Kiyama and M. Itoh, Phys. Rev. Lett. 91, 167202 (2003). 
P. Lunkenheimer, T. Rudolf, J. Hemberger, A. Pimenov, 

S. Tachos, F. Lichtenberg, and A. Loidl, Phys. Rev. B 68, 
245108 (2003). 



12 



11 R. Schmitz, O. Entin-Wohlman, A. Aharony, A. B. Harris, 
and E. Muller-Hartmann, submitted to Phys. Rev. B, cond- 
mat/0407524 (2004). 

12 As compared to Ref. 11, we have slightly modified some 
of the numerical values employed in the calculation of 
the single-bond spin couplings: The effective Ti-O charge- 
tranfer energy used here is A e fi = 6.6 eV instead of the value 
5.5 eV used in Ref. 11. (Hence the spin-exchange couplings 
are multiplied by (5.5/6.6) 2 = 0.694 compared to the ones 
of that reference. This corresponds to a renormalization of 
the units in which the ground-state energy is found, but 
does not affect the minimization of that energy, and con- 
sequently, does not modify the ground-state magnetization 
configuration.) In addition, the spin-exchange couplings are 
computed taking into account only systematic contributions 
in powers of the spin-orbit coupling. 11 



T. Holstein and H. Primakoff, Phys. Rev. 58, 1098 (1940). 
D. Reznik, P. Bourges, H. F. Fong, L. P. Regnault, J. Bossy, 
C. Vettier, D.L. Milius, LA. Aksay, and B. Keimer, Phys. 
Rev. B 53, R14741 (1996). 

S.M. Hayden, G. Aeppli, T. G. Perring, H. A. Mook, and F. 

Dogan, Phys. Rev. B 54, R6905 (1996). 

S. Pailhes, Y. Sidis, P. Bourges, C. Ulrich, V. Hinkov, L. P. 

Regnault, A. Ivanov, B. Liang, C.T. Lin, C. Bernhard, and 

B. Keimer, Phys. Rev. Lett. 91, 237002 (2003). 

B. Keimer, private communication. 

T.A. Kaplan, Z. Phys. B 49, 313 (1983). 

L. Shekhtman, O. Entin-Wohlman, and A. Aharony, Phys. 

Rev. Lett. 69, 836 (1992). 

O. Entin-Wohlman,, A. Aharony, and L. Shekhtman, Phys. 
Rev. B 50, 3068 (1994). 



13 



